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Summary 

The Goddard Cumulus Ensemble model and other models use a technique for simulating 
clouds that includes three types of ice particles: very small cloud ice, larger snow crystals, and 
more dense graupel or hail. The cloud ice particles are important because they produce the large 
ice clouds that are blown off by thunderstorms at high levels in the atmosphere. These cloud ice 
particles grow by the transfer of water vapor onto the particles when the relative humidity is very 
high. As the ice particles grow larger, they eventually become snow. In previous versions of the 
cloud model, this growth process was incorrect because it did not depend on the relative 
humidity of the air. Two corrections for this process are presented that provide the needed 
relative humidity dependence. Computer simulations of thunderstorm lines shows a large impact 
of the corrections. 

Another aspect of the cloud model that is examined is an adjustment that changes excess 
water vapor (i.e., that portion above the amount that the air can hold) into cloud water or cloud 
ice. When the excess water vapor is allowed to form cloud ice, the model produces possible 
errors in the growth of cloud water, ice and snow. It is shown that a significant portion of the 
cloud ice produced by the model within strongly rising air in thunderstorms is the result of this 
adjustment rather than from more realistic cloud processes. 
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ABSTRACT 


Two aspects of the cloud ice parameterization in the Goddard Cumulus Ensemble Model 
cloud physics parameterization are examined: the conversion of cloud ice to snow by depositional 
growth, designated PSFI, and the saturation adjustment scheme. The original formulation of PSFI 
is shown to produce excessive conversion of cloud ice to snow because of an implicit assumption 
that the relative humidity is 100% with respect to water even though the air may actually be quite 
less humid. Two possible corrections to this problem are proposed, the first involving application 
of a relative humidity dependent correction factor to the original formulation of PSFI, and the 
second involving a new formulation of PSFI based on the equation for depositional growth of 
cloud ice. The sensitivity of these formulations of PSFI to the assumed masses of the ice particles 
is examined. Possible problems associated with using a saturation adjustment scheme for cloud ice 
are discussed and simulations of a squall line with and without application of the adjustment 
scheme for ice are compared. 
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1. Introduction 


The Goddard Cumulus Ensemble model (GCE, Tao and Simpson 1993; Tao and Soong 
1986; Tao et al. 1993, 1996) cloud microphysics scheme is patterned after the schemes of Hsie et 
al. (1980; hereafter HFO) and Lin et al. (1983; hereafter LFO), but includes some modifications 
for cloud ice based upon the parameterizations of Rutledge and Hobbs (1983; hereafter RH). The 
parameterization includes three phases of ice: small cloud ice particles, larger snow crystals, and 
more dense graupel or hail particles. The cloud ice particles are important because they contribute 
significantly to the depth, width, and optical thickness of the anvil clouds generated by convection, 
which can significantly impact the transfer of solar and terrestrial radiation in the troposphere. 
Important process in the budget for cloud ice are the initiation of new cloud ice and the conversion 
of cloud ice to snow as the ice crystals grow to the size of snow particles by deposition. 
Examination of some components of the microphysics scheme reveals some potential problems or 
inconsistencies in the parameterization of the cloud ice processes. In this study, we discuss these 
problems and suggest some solutions. 

One of the problems in the LFO and GCE ice schemes is the parameterization of the 
conversion of cloud ice to snow as ice grows by vapor deposition in the presence of cloud water, 
usually referred to as the Bergeron process and designated PSFI. Krueger et al. (1995) pointed 
out that PSFI converts ice to snow even when there is no cloud water. In fact, it converts ice to 
snow independent of the supersaturation with respect to ice and, as a result, produces excessive 
conversion of ice to snow. Since the cloud ice is assumed not to fall while snow does fall, Krueger 
et al. suggested that the excessive conversion of ice to snow via PSFI acts as a crude 
parameterization of cloud ice fallout. This note describes an alternative formulation of PSFI, 
based on the equation for depositional growth in RH, that is dependent on the supersaturation 
with respect to ice. 
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A second potential problem is related to the saturation adjustment scheme (Tao et al. 
1989), which removes supersaturations with respect to water or ice (or both) depending on the 
temperature and on the relative amounts of cloud water and ice. Use of such a saturation 
adjustment scheme can produce cloud water at relative humidities less than saturation and lead to 
possible errors in the growth rates of ice and snow. We will compare simulations with and 
without the saturation adjustment scheme to check for physical consistency of the results. 

2. Numerical Model 

The discussion will be supported with numerical simulations of the 10-11 June 1985 PRE- 
STORM 1 squall line (Johnson and Hamilton 1988; Rutledge et al. 1988; Tao et al. 1993) using the 
two-dimensional version of the GCE model. The model equations are anelastic and the cloud 
microphysics include a parameterized Kessler-type two category liquid, three-category ice (cloud 
ice, snow, graupel/hail) scheme following LFO and Rutledge and Hobbs (1983, 1984). Further 
details about the model can be found in Tao et al. (1993, 1996). 

A stretched vertical grid with 33 grid points and grid spacings from 150 m at low levels to 
1000 m near model top (19 km) was used in order to maximize resolution in the lowest levels of 
the model. The horizontal grid consisted of 1026 grid points, the central 871 of which comprised 
the fine-grid area with a constant 1-km resolution. Outside of this region, the grid spacing was 
horizontally stretched with a ratio of 1.05:1 between adjacent grid points. Open lateral boundary 
conditions were used (Klemp and Wilhelmson 1978). A 5-km deep Rayleigh relaxation 
(absorbing) layer was used at the top of the model. Forward time differencing and a positive- 
definite advection scheme with a nonoscillatory option (Smolarkiewicz and Grabowski 1990) are 
used for all scalar variables (potential temperature, vapor mixing ratio, and all hydrometeor 
categories). A fourth-order accurate advection scheme and leapfrog time integration are used for the 
velocity components. The calculations used a time step of 6 s. 
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The model basic state is derived from the 2330 UTC 1 1 June 1985 sounding taken at Pratt, 
Kansas. This sounding is characterized by convective available potential energy of 2300 J kg '. 
The wind profile has been modified at upper levels to reduce the shear magnitude (see Tao et al. 
1993). Convection is initiated by applying a cooling rate of up to 36 K h' 1 over an area 65 km wide 
and 2.4 km deep for the first 600 s of the simulation. 

Along with changes to the ice microphysical code related to PSFI and the saturation 
adjustment scheme, a few other minor modifications (summarized in Appendix A) have been made 
that generally produce only minor changes in the results. Table 1 summarizes the simulations 
included in this study. Each case is designated by the particular formulation of PSFI (of which 
there are three), the cloud-ice masses used in the parameterization (either values from HFO or from 
RH), and, in some cases, by other letters. The three formulations of PSFI and the definitions and 
values of the ice masses are provided in the next section. For consistency with the equations in 
Koenig (1971) and LFO and the coding of the numerical model, all equations and variables are 
written in cgs units. 

3. Problems in the cloud ice parameterization 

a. Conversion of cloud ice to snow 

The basis for the parameterization of the Bergeron growth process in the ice 
microphysics scheme is Koenig’s (1971) equation for the depositional growth of ice crystals at 
1 00% relative humidity with respect to water. 


dm 

= a.m 

dt 1 


a 2 


( 1 ) 


1 Preliminary Region Experiment for the Stormscale Operational and Research Meteorology Program (Cunning 
1986 ). 
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where m is the crystal mass in grams and a, and a 2 are temperature dependent parameters [see 
Table 4 of Koenig (1971)]. HFO and LFO expressed the conversion of cloud ice to snow through 
vapor deposition as 

PSFI = -^-, (2) 

At, 

where At, is the time required for an ice particle to grow from some initial size m 0 to some 
specified size considered representative of a small snow particle, m s . According to (1), A /, can 
be written as 


At, 


J a , m 2 (1 - a 2 ) 1 J 


m 0 


( 3 ) 


and PSFI can be written as 


PSFI, = 


a,{\-a 2 )q, 


\m 




- 




o 


( 4 ) 


We will refer to (4) as the “original HFO” formulation of PSFI and designate it with the subscript 
1 . 

Following from Koenig’s (1971) formulation for a ] and a 2 , PSFI [ depends only on 
temperature and the specified ice crystal masses m 0 and m s , values for which can be obtained 
from HFO or from other mass-diameter relationships. As currently coded in the models, PSFI t is 
independent of the supersaturation with respect to ice (i.e., 100% relative humidity with respect 
to water is implicit) and overestimates the conversion rate for regions with relative humidities 


4 



less than 100% with respect to water. Furthermore, PSFIi can be active even when the relative 
humidity is less than or equal to the ice saturation value. 

A simple correction to the original HFO formulation of PSFI is obtained by multiplying 
(4) by an empirically derived relative humidity factor (Hindman and Johnson 1972; Reisner et al. 
1998) 


PSFI 2 


(gy-g,.-) a M~ a i) ( i, 

(g w ~ g.v<) - m ^~ a ^ j 


gv > g„- 


( 5 ) 


where q v , q si , and q sw are the vapor mixing ratio, and the saturation vapor mixing ratios over ice 
and water, respectively. Equation (5), which will be referred to as the “modified HFO” 
formulation of PSFI and designated PSFI 2 , allows PSFI 2 to be equivalent to Eq. (4) when the 
relative humidity is 100% and to decrease to zero as q v approaches the ice saturation value. 

Another alternative formulation of PSFI can be derived from the ice depositional growth 
rate equation of RH, 

dm _ 4D.(5, -1) 

dt A" + B" K 

where D i is the mean diameter of the ice crystals, 5, = qjq s , is the supersaturation with respect 
to ice, and A" and B" are essentially temperature dependent parameters given below (A 16) in 
RH. RH assumed hexagonal plate-like crystals and used the mass-diameter relationship 
£) =51.5m l/2 , where m is the average crystal mass in grams and D t is the crystal diameter in 
centimeters. Substituting m for D t in (6) gives 
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( 7 ) 


dm = 206.2(5, - 1 )m l/2 
* “ (A" + B") 


Equation (7) can be expressed in a form similar to (1) by setting «, = 206.2(5, - 1)/( A" + B") and 
a 2 = 0.5 . Equation (4) then gives 


psfi 3 = 


o 

[<-< r 


( 8 ) 


Equation (8) will be referred to as the RH formulation of PSFI and will be designated PSFI3. In 
this formulation, a ] is both temperature and relative humidity dependent and approaches zero as 
the vapor mixing ratio approaches the ice saturation value. Unlike (4)-(5), PSFI 3 is consistent 
with the cloud ice depositional growth rate (6) used in the model. The values of a\ in (8) can be 
compared to the values from Koenig (1971) if we consider relative humidities with respect to 
water of 100% (Figure 1). Koenig’s values reach a maximum value near 258 K and show large 
variations with temperature because of assumptions of different particle habits for different 
temperature regimes. The values of a } from PSFI3 are also maximum near 258 K, but show a 
smoother and much smaller variation with temperature because of the assumption of a single 
particle habit. 

Ice crystal masses from the RH formulation can be substantially different from the 
masses used by HFO for the same assumed size. HFO set m 0 and m s to the masses 

corresponding to 40 micron (m 0 =2.46x1 0 -7 g) and 50 micron (m < =4.80xl0~ 7 g) radius ice 

crystals. Krueger et al. (1995) suggested that 100 microns was a more realistic radius for small 
snow crystals, which gives /w s =3.84xl0~ 6 g. If we use the mass-diameter relationship of RH, 
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then m 0 =2.41xlCT 8 g for a 40 micron radius particle, m 1 =3.76xl0 -8 g for a 50 micron particle, 

and /w s =1.51xl0~ 7 g for a 100 micron particle, all of which are about an order of magnitude 

smaller than the values from HFO. The sensitivity of the numerical results to the varying 
estimates of m Q and m r will be examined later. The designation of HFO or RH in the case name 

indicates which mass values are used in that case. 

Figure 2 shows magnitudes of PSFI based on the three formulations above: the original 
HFO formulation, PSFI | (4), the modified HFO formulation, PSFI 2 (5), and the RH formulation, 
PSFI 3 (8). The calculations assume m 0 and m s corresponding to the masses of 40 and 100- 
micron radius particles using the RH mass-diameter relationship. The cloud ice mixing ratio is 
calculated assuming 40-micron sized particles in concentrations that vary with temperature 
following Fletcher (1962) (Eq. A2), q t =m 0 njp, where p is the air density. Figure 2a clearly 

shows the temperature dependence of the original HFO formulation and the lack of dependence 
upon relative humidity. The rates remain the same even when the air is subsaturated with respect 
to ice. The rates increase with decreasing temperature (increasing height) as a result of the 
increasing number concentration of cloud ice specified by the Fletcher equation. Figure 3 shows 
the snow and cloud ice distributions for case PSFfiRH averaged over hours 5-6 of the simulation 
(using output saved every 1 0 min). The snow field (Fig. 3a) shows an anvil slightly wider than 
150 km, with large mixing ratios contained in the leading convective cells and smaller, more 
horizontally uniform values in the trailing stratiform region. The cloud ice distribution (Fig. 3b) is 
similarly characterized by large values within convective cells, but values in the trailing stratiform 
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region are quite small, a result that will be shown to be caused by excessive conversion of cloud 
ice to snow by PSFI[. 

When the empirically derived relative humidity correction is applied (Fig. 2b, PSFI 2 ), the 
conversion rates for a given temperature are forced to decrease from a maximum value at 100% 
relative humidity with respect to water to zero at the ice saturation value. The RH formulation 
(Fig. 2c, PSFI 3 ) produces a qualitatively similar distribution to the modified HFO formulation, 
but the magnitudes are larger because, at these colder temperatures, the parameter a\ for PSFI3 is 
larger than the Koenig value. Figure 2d shows the difference between the RH formulation and the 
original HFO formulation (PSFI r PSFI 3 ). For ice supersaturations exceeding 10%, the RH 
formulation produces larger growth rates while for small ice supersaturations (<10%) and 
subsaturated conditions, the original HFO formulation produces larger growth rates, particularly 
at colder temperatures. 

The above results imply that PSFI 1 will cause excessive transfer of ice to snow in the 
upper portion of clouds that exist in conditions of near saturation with respect to ice. This fact is 
demonstrated clearly with the numerical results shown in Fig. 4, which depict the snow and cloud 
ice fields for PSFI3 (case PSFI3RH). The snow field (Fig. 4a) shows somewhat reduced mixing 
ratios in the convective region compared to the case using PSFl! (Fig. 3a) while the cloud ice field 
(Fig. 4b) shows substantially greater mixing ratios within most of the anvil, including the trailing 
stratiform region. In particular, PSFI 3 leads to a broadening of the anvil cloud and an increase in 
cloud-top height of about 2-4 km. As in case PSFI[RH, case PSFI3RH results in relative 
humidities near the ice saturation value (Fig. 10b) above ~8.5 km. In case PSF^RH, cloud ice was 
transferred to snow even though depositional growth was negligible, while in case PSFI3RH, no 
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transfer of ice to snow occurred so that cloud ice accumulated in the anvil. If the accumulation of 


ice in the anvil becomes excessive (as determined by either objective or subjective measures), then 
the accumulation of cloud ice can be controlled explicitly through inclusion of cloud ice fallout as 
opposed to having PSFI cause excessive transfer of ice to snow and subsequent fallout of the 
snow. 

To further demonstrate the sensitivity of the cloud ice fields to each formulation of PSFI 
and to the particle masses m 0 and m s . Fig. 5 shows contoured frequency by altitude diagrams 
(CFADs, Yuter and Houze 1995) for cloud ice. The CFAD is a statistical analysis approach 
developed by Yuter and Houze (1995) that succinctly summarizes the frequency distribution of a 
variable in a given volume in a single contour plot. The ordinate of the plot is height, the abscissa 
is the value of the parameter whose distribution is being plotted, and the contours indicate the 
frequency of occurrence per parameter unit per unit height. The CFADs of cloud ice depict the 
frequency of cloud ice amounts within the horizontal domain between 300 to 500 km (e.g., Fig. 
3b) and between hours 4-6 of the simulations. 

Figure 5a shows the CFAD of cloud ice for case PSFIiRH (corresponding to Fig. 3b). 
Mixing ratios greater than 0.2 g kg -1 are seldom produced and most values are less than 0.1 g kg -1 . 
In contrast, the CFAD for case PSFI 3 RH (Fig. 5b) shows greater frequency of cloud ice mixing 
ratios exceeding 0.2 g kg -1 and significantly higher frequency of cloud ice above 10 km. For the 
same mixing ratio amount, the maximum frequency is found at a higher altitude in case PSFI 3 RH 
than in case PSFI|RH. Using the modified HFO formulation (PSFI 2 , Fig. 5c) yields results similar 
to the RH formulation (Fig. 5b). Apparently, the difference in magnitude between PSFI 2 and 
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PSFI 3 (Figs. 2b, c) produces only small changes to the cloud ice fields since the relative humidity 
tends to be near the ice saturation value. 

Previously, it was mentioned that the ice particle masses calculated from the formula of 
RH were about an order of magnitude less than the values assumed in HFO. When the HFO 
particle masses are used instead of the RH values, the factor is about a factor of 

7 larger so that the magnitude of PSFI is reduced by the same factor. Figure 5d shows that when 
the HFO particle masses are used in PSFI i (case PSFIiHFO), the reduced values of PSFI 
contribute to increased ice mixing ratios. When the HFO particle masses are used in PSFI 3 (case 
PSFI3HFO, Fig. 5e), the main differences from case PSFI3RH are seen only between altitudes of 
5-8 km, with little difference apparent above 8 km. The relative lack of sensitivity of PSFI 3 to the 
change in particle masses compared to PSFI! is due to the fact that, at altitudes above 8-10 km, 
the vapor mixing ratios are near or only slightly above the ice saturation value so that PSFI 3 is 
generally small. Since PSFI 1 is independent of supersaturation, a greater reduction of the 
conversion of ice to snow is achieved. 

Another problem inherent in the dependence of PSFI on the assumed masses m s and m 0 
is that as the difference between m s and m 0 becomes small, the conversion rate can become very 
rapid and lead to excessive depletion of cloud ice. If, on the other hand, the difference is large, 
then the conversion rate may become effectively negligible. Figure 6 shows the sensitivity of the 
different formulations of PSFI in Eqs. (4, 5, 8) to the size of the larger threshold particle size, i.e., 
the radius of particles with mass m s , with m s being calculated from the mass-diameter 
relationship of RH. The calculations assume 90% relative humidity with respect to water (120% 
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with respect to ice), air temperature of 243 K, and m 0 equal to the mass of a 40 micron radius 
particle. Significant sensitivity to the radius of the larger particle is seen between 50 and 1 00 
microns. 

Figure 5f (case PSFI 3 RH50) shows that when m s is reduced from that associated with a 
100 micron radius particle to a 50 micron particle, the increase in PSFI 3 significantly reduces the 
cloud ice mixing ratios (compare to Fig. 5b). Somewhat less sensitivity to m s is seen when the 
modified HFO formulation, PSFI 2 , is used (Fig. 5g) because of the smaller exponent, e.g., 
a 2 = 0.436 for T<243 K compared to a 2 = 0.5 for PSFI 3 . Reduced sensitivity is also found when 
the HFO values of m s and m 0 are used (not shown) since these mass values are about an order of 
magnitude larger than in those from the RH mass-diameter formula. Hence, significant sensitivity 
of PSFI to m s is found primarily in the RH formulation (PSFI 3 ) when the RH formula for ice 
mass versus diameter is used. Figure 6 suggests that using m s equal to the mass of a 100 micron 
radius particle brings PSFI 3 into a range where it is less sensitive to small changes in m s . 

Krueger et al. (1995) suggested that the excessive conversion of ice to snow caused by 
PSFI] acts like a crude fallout of cloud ice since ice is assumed not to fall while snow does. This 
suggestion can be tested by including a small fall speed for cloud ice in the model. The 
formulation of PSFI assumes that the cloud ice particles are no larger than about 200 microns in 
diameter. Terminal velocities calculated by Heymsfield (1972) suggest maximum fall speeds of 
20-30 cm s" 1 for particles with lengths of 200 microns. Because of the uncertainty in fall 
velocities and their dependence on crystal habit, we simplify matters by using a constant fall 
velocity of 20 cm s' 1 in the model. Figure 7 shows the snow and cloud ice fields for case 
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PSFI3RHVI, which includes cloud ice fallout. The 20 cm s' 1 fall speed has only a small effect on 
the snow and cloud ice distributions, the main effect being a lowering of the cloud top by about 1 
km (cf. Figs. 4 and 7). Comparison of Fig. 7 to Fig. 3 also suggests that the excessive conversion 
of ice to snow by PSFIi is not equivalent to cloud ice fallout. 

b. The saturation adjustment scheme 

The saturation adjustment scheme in the GCE code is designed to account for 
condensation at warm temperatures (T> 273 K), deposition at temperatures less than 233 K, and 
a linear mix of the two at intermediate temperatures. Three problems exist with such a scheme. 
First, for 233 K < T < 273 K, cloud water can be produced even when the relative humidity with 
respect to water is less than 100%. Second, for T < 273 K, the scheme produces deposition in 
addition to the explicit calculation of deposition in the term PIDEP (see RH), if supersaturation 
with respect to ice is still present after the calculation of PIDEP and ice initiation, PINT 
(Appendix A). In particular, the effect is to produce more cloud ice at middle levels (243 K < T < 
273 K) than is parameterized by PIDEP. Finally, since the saturation adjustment scheme is 
computed prior to the depositional growth of snow and graupel, it reduces the depositional 
growth of the precipitation categories. This problem can be particularly troublesome in stratiform 
regions where significant growth of snow by deposition is known to occur (Rutledge and Houze 
1987; Braun and Houze 1994; Houze 1997). 

Figure 8 shows the snow and cloud ice distributions for case PSFI 3 RHNS, which does not 
include the saturation adjustment with respect to ice. Snow mixing ratios (Fig. 8a) are comparable 
to the case with the saturation adjustment for ice (PSFI 3 RH, Fig. 4a), but with slightly larger 
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snow mixing ratios in the trailing stratiform region. The cloud ice mixing ratios (Fig. 8b) are 
significantly reduced in the convective region and slightly reduced in the stratiform region. In 
particular, very little cloud ice is present below 8 km without the saturation adjustment for ice. 
These changes are quite evident in the CFAD for cloud ice (Fig. 5h), which shows mixing ratios 
rarely exceeding 0.1 g kg -1 above 8 km and 0.01 g kg' 1 below. It is evident that a significant 
portion of the cloud ice is produced by the saturation adjustment scheme and not by the 
parameterized ice microphysical processes (PINT, PIDEP). 

As mentioned above, one problem with the saturation adjustment scheme, when applied 
to cloud ice, is that it can produce cloud water when the relative humidity is less than 100% with 
respect to water. Figure 9 shows cloud water mixing ratios for cases PSFI 3 RFI and PSFI3RHNS. 
Figure 9a indicates that cloud water amounts for case PSFI3RH occasionally exceed 0.1 g kg -1 in 
regions that have relative humidities less than 98%, with these near saturated regions extending 
up to about 5.5 km. Figure 9b shows a qualitatively similar cloud water distribution in case 
PSFI3RHNS, but in this case, cloud water is virtually always contained within regions having 
relative humidities exceeding 98%, with these near saturated regions extending up to about 8 km. 

The relative humidities with respect to ice in cases PSFI3RH and PSFI3RHNS are quite 
different. For case PSFI3RH (Fig. 10a), the relative humidity field shows very uniform values 
near the ice saturation threshold within the anvil cloud above 9-10 km and relatively small 
supersaturations below 9 km. The relative humidity is very uniform at levels where T < 233 K 
since the adjustment scheme brings the vapor mixing ratio to its saturation value, which is simply 
a function of temperature and has little small-scale variability within the upper portion of the 
cloud. Tendrils of cloud ice and high humidity hang down from the upper cloud canopy and 
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result from the saturation adjustment rather than from the explicit calculation of ice deposition in 
the term PIDEP. Since the supersaturations with respect to ice are small, the conversion rate 
PSFI 3 is also small. In contrast, in case PSFI 3 RHNS (Fig. 10b), the relative humidity field shows 
more variability at upper levels and contains very high supersaturations between 6-8.5 km. There 
is no longer a constraint that the relative humidity be near the ice saturation value, but above 9 
km, the relative humidity is still generally very close to this value. The tendrils of cloud ice are 
not apparent in this case since the ice deposition term PIDEP, being a strong function of 
temperature because of the Fletcher equation (Eq. Al), is unable to produce substantial amounts 
of ice at temperatures warmer than about 243 K. Furthermore, where the relative humidities with 
respect to ice are large, PSFI 3 rapidly converts the ice to snow. 

The reduction in supersaturation by the saturation adjustment scheme (involving both 
cloud water and ice) before the calculation of snow depositional growth substantially reduces the 
potential growth of the snow field. Turning off the saturation adjustment for ice increases the 
relative humidity with respect to ice so that snow grows more by deposition, but it also reduces 
the amount of cloud ice so that snow grows less by conversion of ice to snow via PSFI and by 
accretion of cloud ice. This latter process, however, is highly dependent upon the collection 
efficiency of cloud ice by snow, the value of which is uncertain and probably overestimated in 
the model [see Eq. (23) and Fig. 4 of Lin et al. (1983)]. Consequently, turning off the saturation 
adjustment for ice may not necessarily lead to significant changes in snow production. 

4. Conclusions 

This note has described some potential problems associated with a parameterization of 
cloud ice processes that is used in many fine-scale cloud models, including the Goddard Cumulus 
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Ensemble (GCE) model. One problem is associated with the rate at which cloud ice is transferred 
to snow as a result of growth by vapor deposition, designated PSFI,. The formulation generally 
used in the parameterization is independent of relative humidity, which causes ice to be converted 
to snow even when the air is subsaturated with respect to ice. Two alternative formulations are 
presented. The first alternative, designated PSFI 2 , simply multiplies the original formula by an 
empirically derived relative-humidity dependent factor so that PSFI 2 diminishes as the relative 
humidity approaches the ice saturation value. The second alternative formulation, designated 
PSFI 3 , is derived directly from the equation for depositional growth of cloud ice used in the model 
(PIDEP). This formulation causes PSFI, to diminish as the relative humidity approaches the ice 
saturation value, but also ensures physical consistency with the parameterization of depositional 
growth of cloud ice. The two alternative formulations produce relatively similar results. Their main 
impacts include an increase of cloud-top height and a substantial increase in the cloud ice mixing 
ratios, particularly at upper levels in the cloud. 

A second issue of concern in the cloud ice parameterization is the use of a saturation 
adjustment technique for ice. While the technique is useful for generating cloud ice at temperatures 
warmer than about 243 K, it has the disadvantages in that it produces cloud water in slightly 
subsaturated air, forces cloudy air at cold temperatures (< 243 K) to have relative humidities at the 
ice saturation value, and prevents large supersaturations with respect to ice, which then limits the 
growth of snow by vapor deposition. Simulations in which the saturation adjustment was applied 
only to cloud water resulted in substantially reduced cloud ice amounts in the convective regions 
and nearly eliminated cloud ice at temperatures warmer than -243 K. These results suggest that the 
saturation adjustment technique was responsible for generating a substantial portion of the cloud 
ice in convective cells and at temperatures where mixed-phase processes are important. A more 
realistic treatment of cloud ice processes should probably include cloud ice concentration as a 
prognostic variable in order to reduce the dependence of the ice scheme on the Fletcher equation for 
the number concentration of ice nuclei and should include improved formulations of the processes 
that produce cloud ice including ice multiplication mechanisms (e.g., Reisner et al. 1998). 
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APPENDIX A 

The number concentration of ice nuclei is specified according to Fletcher (1962) as 


n c =n 0 exp[p(T 0 -T)\, 


(Al) 


where T is temperature, 7" 0 =273 K, and [5 and n 0 are often taken as 0.6 and 10' 8 cm' 5 , 

respectively 2 . One problem with the use of (Al) is that at very cold temperatures ( T < 243 K), n c 
attains unrealistically high values, which impacts the calculation of the initiation of new ice particles 
(PINT) and the depositional growth of existing ice particles (PIDEP). To avoid overestimates of 
n c , two alternative solutions are 1) to limit ( T 0 -T) to some specified temperature deficit, say 31 K, 
or to limit the maximum concentration of ice to some prescribed value such as 1 cm' 5 . The impact 
of these solutions is small when the saturation adjustment scheme is used but can be significant 
when the saturation adjustment scheme is not used (not shown). 

The initiation of cloud ice, PINT, was formulated by Rutledge and Hobbs (1983) as 


2 We have changed /3 to 0.46 and n c to 10” 6 cm” 3 in the simulations in this study because these values gave 

concentrations that were in better agreement with some observed values shown in Fig. 1 of Meyers et al. (1992). 
Tests with the model indicate only minor changes in the cloud ice mixing ratios produced by the changes in /3 and 


16 



PINT = min 


M 0 w t . q v - q si 

pAt ’ At 


(A2) 


where n c is the prescribed number of ice particles, M 0 is the approximate mass of initiated ice 
particles, pis the density of air, and At is the model time step. This production term is assumed to 

produce n c particles with mass M 0 , unless it removes all of the existing supersaturation. At many 
grid points, cloud ice already exists when a supersaturation is produced. Rather than initiating n c 
particles, PINT should initiate n c = n c -n ( particles ( n= pq/M ], M,=average mass of existing ice 
particles), i.e., it should only initiate the number of particles needed to bring the concentration up 
to the prescribed number, n c , possible at the ambient temperature. An alternative formulation of 
(A2) can be written as 


PINT = min 


M o q v - 
pAt ’ At 


(A3) 


Equation (A3) is similar to the formulation used by Reisner et al. (1998). Model results indicate 
that replacing (A2) with (A3) in the microphysical scheme produces only relatively small changes 
in the cloud ice field. 
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Table 1. Summary of the numerical simulations. The column labeled “ice mass” indicates whether 


the HFO or RH values are used for the initial and final sizes of ice crystals, m 0 and m s , the values 


of which are provided in the last two columns. Case PSFI,RHVI differs from PSFI,RH in that it 
includes a fall velocity for cloud ice. Case PSFI,RHNS differs from PSFI,RH in that it does not 


use the saturation adjustment for cloud ice. 


Name 

PSFI 

Formulation 

Ice 

Mass 

™ o 

(g) 

m s 

(g) 

PSFI.RH 

PSFI, 

RH 

2.4lxl0' 8 

1.51xl0' 7 

PSFI 2 RH 

PSFI, 

RH 

2.41xl0" 8 

1.51xl0' 7 

psfi,rh 

psfi; 

RH 

2.41xl0‘ 8 

1.51xl0 7 

PSFI.HFO 

PSFI, 

HFO 

2.46xl0' 7 

3.84x1 O' 6 

PSFI 2 HFO 

psfi 2 

HFO 

2.46xl0' 7 

3.84x10 6 

psfi,hfo 

PSFI, 

HFO 

2.46xl0' 7 

3.84x10 6 

PSFI ? RH50 

psfi 2 

RH 

2.41xl0‘ 8 

3.76x10 8 

PSFI,RH50 

PSFI, 

RH 

2.41xl0' 8 

3.76xl0' 8 

psfi,rhvi 

PSFI, 

RH 

2.41x10 8 

1.51xl0' 7 

PSFI,RHNS 

PSFI, 

RH 

2.41xl0‘ 8 

1.51x10 7 
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FIGURE CAPTIONS 


Figure 1. Curves of a\ as a function of temperature. The solid line shows the values from Koenig 
(1971) while the dashed line shows the values for the new formulation of PSFI 3 . For the latter 
values, a relative humidity of 100% has been assumed in order to be consistent with Koenig. 

Figure 2. Magnitude (in units of 10~ 7 g g _l s _l ) of the various formulations of PSFI: a) the original 
HFO formulation (Eq. 4), b) the modified HFO formulation (Eq. 5), and c) the RH formulation 
(Eq. 8). Temperature is shown along the vertical axis of the plots as a proxy for height. For these 
calculations, the temperature was assumed to vary in the vertical according to the 2330 UTC 
Pratt, Kansas, sounding. Solid lines in (a-c) indicate conversion rates, thick dashed lines 
supersaturation with respect to ice. (d) Difference between PSFI] and PSFI 3 , with positive 
(negative) values indicated by solid (dashed) lines. 

Figure 3. Vertical cross sections of time-averaged (a) snow mixing ratio and (b) cloud ice mixing 
ratio for case PSFI]RH. Cross sections were obtained by averaging fields from hours 5-6 using 
output at 10-min intervals. The contour intervals are 0.1 g kg' 1 starting at 0.01 g kg' 1 for snow 
and 0.025 g kg' 1 starting at 10' 3 g kg' 1 for cloud ice. 

Figure 4. Same as in Fig. 3, but for case PSFI 3 RH. 
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Figure 5. Contoured Frequency by Altitude Diagrams (CFADs) for cloud ice mixing ratio. 
Contours indicate the frequency km' 1 (g kg' 1 )' 1 . Contours are drawn at values of 0.5, 1, 2, 4, 8, 
12, 20, 36, 48, 60, and 72. The 4 km' 1 (g kg" 1 )' 1 contour is highlighted by the dashed line. 

Figure 6. Sensitivity of PSFI to the threshold radius of snow corresponding to the mass m s in 
Eqs. (4, 5, and 8). 

Figure 7. Same as in Fig. 4 (case PSFI3RH), but for case PSFI3RHVI, which includes a fall 
velocity for cloud ice. 

Figure 8. Same as in Fig. 4 (case PSFI3RH), but for case PSFI3RHNS, which does not use the 
saturation adjustment for cloud ice. 

Figure 9. Vertical cross sections at t = 6 h of cloud water (contours) and relative humidity with 
respect to water (shading) for cases (a) PSFI3RH (with the saturation adjustment for ice) and (b) 
PSFI3RHNS (without the saturation adjustment for ice). Shading indicates regions with relative 
humidity exceeding 98%. Contours for cloud water are drawn at values of 0.01, 0.1, 0.5, 1, 2, and 
3 g kg' 1 . 

Figure 10. Vertical cross sections at t= 6 h of cloud ice (contours) and relative humidity with 
respect to ice (shading) for cases (a) PSFI3RH (with the saturation adjustment for ice) and (b) 
PSFI3RHNS (without the saturation adjustment for ice). Shading indicates regions with relative 
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humidity exceeding 100% (light shade), 110% (dark shade), and 120% (white). Contours for 
cloud ice are drawn at values of 0.01, 0.05. 0.1, 0.15, 0.2, and 0.3 g kg' 1 . 
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Figure 1. Curves of a] as a function of temperature. The solid line shows 
the values from Koenig (1971) while the dashed line shows the values for 
the new formulation of PSFI 3 . For the latter values, a relative humidity of 
100% has been assumed in order to be consistent with Koenig. 
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Figure 2. Magnitude (in units of 10" 7 g g' 1 s -1 ) of the various formulations of PSFI: 

(a) the original HFO formulation (Eq. 4), (b) the modified HFO formulation (Eq. 5), 
and (c) the RH formulation (Eq. 8). Temperature is shown along the vertical axis of 
the plots as a proxy for height. For these calculations, the temperature was assumed to 
vary in the vertical according to the 2330 UTC Pratt, Kansas, sounding. Solid lines in 
(a-c) indicate conversion rates, thick dashed lines supersaturation with respect to ice 
(qMsi)- (d) Difference between PSFI, and PSFI 3 , with positive (negative) values indi- 
cated by solid (dashed) lines. 
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Figure 3. Vertical cross sections of time- 
averaged (a) snow mixing ratio and (b) 
cloud ice mixing ratio for case PSF^RH. 
Cross sections were obtained by averaging 
fields from hours 5-6 using output at 10- 
min intervals. The contour intervals are 0. 1 
g kg ' 1 starting at 0.01 g kg " 1 for snow and 
0.025 g kg ' 1 starting at 0.001 g kg ' 1 for 
cloud ice. 


Figure 4. Same as in Fig. 3, but for case 
PSFI 3 RH. 
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Figure 5. Contoured Frequency by Altitude Diagrams (CFADs) for cloud ice mixing 
ratio. Contours indicate the frequency km' 1 (g kg" 1 )' 1 . Contours are drawn at values 
of 0.5, 1, 2, 4, 8, 12, 20, 36, 48, 60, and 72. The 4 km' 1 (g kg' 1 )' 1 contour is highlighted 
by the dashed line. 

















Figure 6. Sensitivity of PSFI to the threshold radius of 
snow corresponding to the mass m s in Eqs. (4,5, and 8). 
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Figure 7. Same as in Fig. 4 (case PSFI 3 RH), but 
for case PSFI 3 RHVI, which includes a fall 
velocity for cloud ice. 







Figure 8. Same as Fig. 4 (case PSFI 3 RH), but for 
case PSFI3RHNS, which does not use the saturation 
adjustment for cloud ice. 
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Figure 9. Vertical cross sections at t=6 h of cloud water (contours) and rela- 
tive humidity with respect to water (shading) for cases (a) PSFI 3 RH (with the 

saturation adjustment for ice) and (b) PSFI 3 RHNS (without the saturation 
adjustment for ice). Shading indicates regions with relative humidity exceed- 
ing 98%. Contours for cloud water are drawn at values of 0.01, 0.1, 0.5, 1 , 2, 
and 3 g kg' 1 . 






300 350 400 450 500 

Distance (km) 


Figure 10. Vertical cross sections at t=6 h of cloud ice (contours) and relative 
humidity with respect to ice (shading) for cases (a) PSFI 3 RH (with the satura- 
tion adjustment for ice) and (b) PSFI 3 RHNS (without the saturation adjust- 
ment for ice). Shading indicates regions with relative humidity exceeding 
100% (light shade), 1 10% (dark shade), and 120% (white). Contours for cloud 
ice are drawn at values of 0.01, 0.05, 0.1, 0.15, 0.2, and 0.3 g kg' 1 . 




